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Abstract 

Using Euler's formula for a network of polygons for 2D case (or 
polyhedra for 3D case), we show that the number of dynamic degrees 
of freedom of the electric field equals the number of dynamic degrees 
of freedom of the magnetic field for electrodynamics formulated on a 
lattice. Instrumental to this identity is the use (at least implicitly) 
of a dual lattice and of a (spatial) geometric discretization scheme 
based on discrete differential forms. As a by-product, this analysis also 
unveils a physical interpretation for Euler's formula and a geometric 
interpretation for the Hodge decomposition. 
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1 Introduction 



There has been continual interest in formulating physical theories on a dis- 
crete lattice [1] [2] [3] [4] . This "latticization" process not only ehminates 
the infinities but also provides a computational route to study the non- 
perturbative region, which it is often difficult to handle by purely analytical 
tools. In this letter, we shall discuss aspects of a lattice version of classical 
electrodynamics based on a geometric discretization [5] [6] . 

There are three main approaches to discretize Maxwell equations: finite 
differences [7] [8], finite elements [9] and finite volumes [7]. However, naive 
implementations of these finite methods to discretize Maxwell equations on 
irregular lattices are known to cause problems, such as spurious modes [10] 
and late-time unconditional instabihties [11], that destroy the solutions. This 
is often a consequence of the failure to capture some essential physics of the 
continuum theory. 

By applying some tools of algebraic topology and a discrete analog of 
differential forms, classical electrodynamics can be constructed from first 
principles on a lattice [5]. The purpose of this paper is to show that via 
such geometric discretization, the equality (necessary, but often not trivially 
observed in common discretization schemes) between the number of dynamic 
degrees of freedom {DoFs) of the electric field and the number of dynamic 
DoFs of the magnetic field is a simple consequence of Euler's formula for a 
network of polygons for 2-1-1 case or volume polyhedra in the 3-1-1 case. As 
a by-product, this analysis also unveils a physical interpretation for Euler's 
formula and a geometric interpretation for the Hodge decomposition. 

2 Lattice electromagnetic theory 

For simplicity, wc consider first a TE field in 2+1 dimensions. The extension 
to 3-1-1 dimensions is considered in Section 3.2. 
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Maxwell equations in the Fourier domain are written in terms of differ- 
ential forms as [5] 

dE = iuB 

dB ^ 

dH = -iojD + J 

dD ^ Q (1) 

In the 2+1 case, H is a, 0-form, D , J and E are 1-forms, and B and Q are 
2-forms. The operator d is the exterior derivative. Prom the identity d'^ — 0, 
the electric current density J and the charge density Q satisfy the continuity 
equation 

dJ = iuQ (2) 

Constitutive equations, which include all metric information in this frame- 
work, are written in terms of Hodge star operators (that fix an isomorphism 
between p forms and 2 — p forms in the 2-1-1 case) [5] 



D = ^eE ,B = (3) 

In numerical simulations, because of limited memory, the computational 
domain is taken as a closed domain D, with boundary dQ. In the lattice, the 
boundary dQ is approximated by a set of linked edges dQ (see Fig.l). 

The latticization corresponds to tiling Q with a finite number Np of poly- 
gons S^, m — 1, .., Np, of arbitrary shape(see Fig. 2) 

^ Np 

Q ~ Q = © S„ (4) 

m=l 

We require the tiling to be conformal i.e., two polygons are either connected 
by one single edge or are not connected at all (see Fig.3). These polygons 
should also be oriented, forming a cell-complex (see Fig.4) [5] [12] .We denote 
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such oriented tiling (Fig.4) the primal lattice. Prom the primal lattice, one 
can construct a dual lattice by connecting interior points of each adjacent 
polygon ^. The dual lattice inherits an orientation from the primal lattice. 

Now we consider casting Maxwell equations on a lattice using the natural 
latticization provided by casting differential forms of various degrees p in 
Eq. (1) as dual elements (cochains) to p dimensional geometric constituents 
of the lattice, i.e., cells: nodes, edges and faces [5]. In the primal lattice, 
we associate the electrostatic potential (0-form ) with primal nodes (0- 
cells), the electric field intensity E (1-form) with primal edges (1-cells) and 
the magnetic flux density B (2- form) with primal cells (2-cells). In the dual 
lattice, we associate the magnetic field intensity H (0-form ) with dual nodes 
(0-cells), the electric flux density D (1-form), the electric current density J 
(1-form) with dual edges (1-cells), and the charge density Q (2-form) with 
dual cells (2-cells). This is illustrated in Fig.5. 

The exterior derivative d can be discretized via its adjoint operator, the 
boundary operator, d, by applying the generalized Stoke's theorem on each 
p-cell of the cell-complex 

{Y^\da^)^{dY^\a^) , (5) 

where 7^+^ is a p+1 dimensional cell and G F^{Q), with (Q) being 

the space of differential forms (cochains in the discrete setting) of degree p 

on the domain Q. We denote an ordered sequence of the above pairing of 

cochains with each of the cells by block letters E, B, H, D, J, Q in what 

follows. These are the DoFs of the lattice theory. In terms of these DoFs, 

^Different geometrical constructions can be used to define the dual lattice. We will not 
dwell here into such discussion since our conclusions depend only on topological properties 
(connectivity) of the lattice. 



4 



the lattice analog of Maxwell's equations is written as [5] [13] 



CE = icuM (6) 

= (7) 

CU = -iuB) + S (8) 

SB ^ Q (9) 



In the above, C,C,S,S are the incidence matrices [5], obtained by applying 
the boundary operator d on each p-cell (discrete version of the exterior deriva- 
tive d). Because the exterior derivative d is a purely topological operator 
(metric-free), the incidence matrices represent pure combinatorial relations, 
whose entries assume only {—1,0, 1} values. 

The lattice version of the Hodge isomorphism can be, in general, written 
as follows 

D=K]E, B=[*^]H (10) 

where both [kg] and [k^] are square invertible matrices. We will not discuss 
here how to construct the Hodge matrices [^e] and [*^]. That the metric 
dependent Hodge matrices are the square invertible matrices is the only as- 
sumption needed here. 

In many problems of interest, one is usually interested only in the dynamic 
behavior of the field, which is determined by the dynamic DoFs of lattice 
theory. Here we show that the number of dynamic DoFs of the electric field, 
DoF'^[E), equals to the number of dynamic DoFs of the magnetic flux, 
DoF'^ (B), in the discretization above. Furthermore, from the isomorphisms 
between E and D, and between H and B (from the Hodge maps), this also 
implies 

DoF"^ {E) = DoF'^ (B) = DoF'^ (D) = DoF'^ (H) (11) 

where the superscript d stands for dynamic. Moreover, for 2-1-1 TE field, this 
number equals to the total number of polygons used for tiling the domain Q 
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minus one {Np — 1). 



3 Hodge decomposition 



The Hodge decomposition can be written in general as 



pp (n) = dpp-^ (Q) e 5FP+^ (n) e {^) 



(12) 



where (^) is the finite dimensional space of harmonic forms, and 5 is the 
codifferential operator, Hilbert adjoint of d [14]. Applying (12) to the electric 
field intensity E, we obtain 



where is a 0-form and A is a 2- form. In Eq. (13) dcf) represents the static 
field and SA represents the dynamic field, and x represents the harmonic 
field component. 

3.1 2+1 theory in a contractible domain 

If domain Q is contractible, x is identically zero and the Hodge decomposition 
can be simplified to 



In the present lattice model, the number of DoFs for the zero eigenspace 
{cu = 0) equals the number of internal nodes of the primal lattice [10] [15] 
[16]. This is because the DoFs of the potential 0, which is a 0-form, is 
associated to nodes. 

Now we show the identity (11). Recall the Euler's formula for a general 



E ^ d(f) + 5A + x 



(13) 



E^d(t) + SA 



(14) 
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network of polygons without holes (Fig.2) 



Nv-Ne^1-Nf (15) 

Here, Ny is the number of vertices (nodes), A^^; the number of edges, and 
Np the number of faces (cells). For any 90, it is easy to verify that 

K-N'e = (16) 

where Ny is the number of vertices on the boundary and the number of 
edges on the boundary (the superscript b stands for boundary). Note that 
cochains on dO, are not associated to DoFs, since they are fixed from the 
boundary conditions. Using the Hodge decomposition (13), the number of 
dynamic [uj ^ Q ) DoFs of the electric field, corresponding to 5A, is given 

by 

DoF'^{E) = N'^-N^ 

= NE-Ny (17) 

where the superscript in stands for internal. Since E is given along the 
boundary, then, for uj ^ 0, J^B is fixed by 

iu I I E (18) 
Jo, Jan 

This corresponds to one constraint on B. Subtracting one degree of freedom 
from the constraint (18), the number dynamic DoFs of the magnetic fiux B 
is 

DoF'^ (B) ^Nf-1 (19) 
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Prom Euler's formula (15), we then have the identity 

DoF'^ (E) = DoF'^ (B) (20) 

Furthermore, thanks to the Hodge isomorphism, the identity (11) follows 
directly 

3.2 3+1 theory in a contract ible domain 

The source free Maxwell equations in 3+1 dimensions read as 

dE = iuB, (21) 

dB = 0, (22) 

dH = -iLvD, (23) 

dD = (24) 

where now H and E are 1-forms, and D and B are 2-forms. The spatial 
domain Q is again (approximately) tiled by a set of polyhedra Q and the 
boundary is by a polyhedron dCl . Using Euler's formula for fi, we have 

Ny-Nj,^l-Np + Np (25) 
and Euler's formula for the boundary polyhedron dO, 

N^-N^E^2-N^p (26) 

where Np is now the number of polyhedra. Combining Eq. (25) and (26), 
we obtain 

{Ne - iV|) - {Ny - K) = {Np - - (iVp - 1) (27) 
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Using the Hodge decomposition (14), the number of dynamic DoFs of the 
electric field (corresponding to 5 A ) is 

= {Ne - Nl;) - {Nv - N'y) (28) 

Each polyhedron produces one constraint for the magnetic flux B from 
Eq.(22). Furthermore, this set of constraints span the condition at the bound- 
ary dn. The total number of the constrains for B is therefore {Np — 1) . 
Consequently, the number of DoFs for the magnetic flux B is 

DoF'^{B) = Np-{Np-l) 

= {Np - N'p) - {Np - 1) (29) 

Identity (11) then follows from Eq. (27), (28) and (29). 

3.3 2+1 theory in a non-contractible domain 

Now consider a non-contractible two-dimensional domain Q with a finite 
number g of holes (genus). This is illustrated in Fig. 6 for 5- = 1. Along the 
boundary of each hole, the electric field E is constrained by 

J E^M (30) 

where the magnetic current density ^ M (passing through the hole) is a 

known quantity. The equation (30) accounts for the possible existence of 

^In physical terms, the magnetic current density M is identified with the "displace- 
ment magnetic current density" iujB, which is given for some cases. In some other cases, 
M comes also from the equivalent magnetic current density by the surface equivalence 
theorem [17]. It should be emphasized that the equivalent magnetic current results from 
an impressed electric field E, not from the movement of any "magnetic charge". 
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the harmonic forms x on In particular, the number of holes g is equal 
to the dimension of the space of harmonic forms x ^^id gives the number 
of independent constraint equations (30). Subtracting g from Eq. (17), the 
number of dynamic DoFs of the electric field in this case becomes 

DoF'^{E) = N'^-N^-g 

= N^-Nv-g (31) 

whereas the number of DoFs of the mag netic flux DoF'^ (B) remains Np — 1 . 
Since Euler's formula for a network of polygons with g holes is 

Nv-NE = {l-g)- Nf (32) 

we have that from Eq. (19), (31) and (32), the identity (11) is again satisfied. 

3.4 Euler's formula and Hodge decomposition 

From the above considerations, we can trace the following correspondence 
in the 2+1 case 

Ne ^ Nv + {Nf-1) + g (33) 

II II 
E ^ d(j) + SA + X 

The number of edges Ne corresponds to the dimension of the space of (dis- 
crete) electric field intensity E (1-forms), which is the sum of the number of 
nodes Ny (dimension of the space of discrete 0-forms 0) , the number of faces 
(Np — 1) (dimension of the space of discrete 2-form A) and the number of 
holes g (dimension of the space of harmonic form x) ■ These correspondences 
attach a physical meaning to Euler's formula and a geometric interpretation 
to the Hodge decomposition. We note that the identity (33) can be viewed 
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as 

N'i = Np + {Nf-1) + g (34) 

II II 
E ^ dcf) + 5A + X 

since only the internal edges and nodes describe the degrees of freedom. We 
can simply drop the superscript in because the identity (16). For the 3+1 
similar correspondence could also be drawn. 

4 Concluding remarks 

Based on a geometric discretization^, we have shown that Euler's formula 
matches the algebraic properties of the discrete Helmholtz decomposition in 
an exact way. Furthermore, we have showed that the number of dynamic 
DoFs for the electric field equals the number of dynamic DoFs for the 
magnetic field on such lattices'^. 

Regarding the time discretization, we also remark that simplectic inte- 
grators, originally developed for Hamiltonian systems [22], can provide a 
time discretization that respects the simplectic structure [23]. However, it is 
not trivial to formulate the Maxwell's equations on a lattice as the canonical 
equations of the Hamiltonian, because electrodynamics can be thought of as 
a constrained dynamic system. A Hamiltonian requires that the canonical 
pair {E, H) should have the same number of degrees of freedom. Identity 
(11) suggests that it is indeed possible to formulate electrodynamics on a 
lattice as the canonical equations of the Hamiltonian. 

^One key feature of this scheme is the use of a dual lattice and of a geometric discretiza- 
tion scheme based on differential forms, also proposed in different contexts in [18] [19] [20]. 

^For the case of high order 1-forms [21], the DoFs of 1-forms could associate with the 
faces and volumes. However, the dimension of the range space of 1-forms (e.g. E) equals 
the dimension of 2-forms (e.g. B) if such 2-forms have zero range space due to dB = 0, so 
the identity (11) still holds. 
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Figure captions 



Fig.l. The curved boundary dfl is approximated by a set of linked edges 

on. 

Fig. 2. Tiling the computation region with polygons. 

Fig. 3 (a) conformal tiling (cell complex); (b) non-conformal tiling. 

Fig. 4. Oriented polygons. 

Fig. 5. Solid lines represent the primal lattice. Primal nodes (vertices) 
are paired with (e.g., node 1), primal edges with E (e.g., edge 15) and 
primal cells with B (e.g., cell 12345). Dashed lines represent the dual lattice. 
Dual nodes are paired with H (e.g., node 4'), dual edges with {D, J) (e.g., 
edge 3'4') and dual cells with Q (e.g., cell 1'2'3'4'5'). 

Fig. 6. 2+1 theory in a non-contractible domain (network of polygons 
with a hole, illustrated by a triangle 123). 



15 



This figure "Fig.l.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/hep-lat/0408005v2 



This figure "Fig. 2.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/hep-lat/0408005v2 



This figure "Fig. 3.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/hep-lat/0408005v2 



This figure "Fig.4.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/hep-lat/0408005v2 



This figure "Fig. 5.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/hep-lat/0408005v2 



This figure "Fig. 6.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/hep-lat/0408005v2 



